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Abstract 

kAn  experimental  investigation  was  conducted  to  determine  the 
validity  of  an  algorithm  developed  by  Whaley  (Ref  14)  to  approximate 
the  natural  frequencies  of  a-  complex  structure  under  arbitrary  mass 
loading  conditions,  when  only  the  unloaded  natural  frequency  and 
mode  shape  data  is  known.  The  chosen  test  specimen  was  a  curved, 
rib-stiffened  panel  from  a  C-141  Starlifter,  aircraft.  The 
panel  was  suspended  from  the  ceiling  by  bungy  coir  and  tested 
in  an  unloaded  configuration  and  nine  separate  mass  loaded  config¬ 
urations.  Then  using  only  unloaded  data  the  generalized  mass  and 
generalized  stiffness  for  each  mass  loaded  configuration  were  computed, 
and  the  natural  frequencies  for  each  configuration  were  computationally 
predicted  using  the  aforementioned  algorithm.  The  theoretical  and 
experimental  results  were  then  compared  to  determine  the  amount  of 
error  incurred  in  the  approximation  technique.  The  theory  of  how 
to  ultimately  determine  the  overall  forced  response  of  the  specimen 
was  discussed  and  an  error  model  was  developed  to  enable  an  examination 
of  the  reliability  of  the  algorithm  in  predicting  forced  response. 
Recommendations  concerning  future  test  procedures,  areas  requiring 

further  study,  and  the  use  of  the  algorithm  were  made. i 

"  \ 


\ 


ix 


THE  PREDICTION  OF  MASS  LOADED  NATURAL  FREQUENCIES 
AND  FORCED  RESPONSE  OF  COMPLEX,  RIB-STIFFENED  STRUCTURES 


I  Introduction 


Background 

The  thrust  of  this  study  is  aimed  at  the  ever-increasing  need  for 
better  performance  in  the  field  of  vibration  control  methods  for  airborne 
electro-optical  systems.  In  order  to  yield  acceptable  performance  levels 
in  such  a  scheme,  the  designer  must  be  provided  with  an  accurate  descrip¬ 
tion  of  the  vibration  characteristics  of  the  airframe.  Also,  these 
vibration  characteristics  may  change  markedly  due  to  the  additional 
mass  of  the  electro-optical  system.  This,  due  to  the  sensitivity  of 
the  electro-optical  system  to  airframe  vibration  disturbances,  may  cause 
serious  degradation  of  system  accuracy.  Therefore,  one  would  also  like 
to  be  able  to  predict  how  the  vibration  characteristics  are  altered  due 
to  arbitrary  mass  loading  conditions.  In  an  effort  to  make  this  possible, 
several  approximate  methods  for  prediction  of  the  modified  characteristics 
have  been  developed.  Some  of  these  are  discussed  later  in  this  chapter. 
The  method  considered  in  this  study  was  developed  by  Whaley  (Ref  14) . 

The  algorithm  he  developed  will  be  derived  later  in  this  study.  It 
allows  a  designer,  given  the  unloaded  natural  frequencies  and  mode 
shapes  of  a  particular  structure,  to  predict  the  change  in  the  natural 
frequencies  due  to  any  arbitrary  mass  loading  conditions.  Tests  have 
been  conducted  on  a  cantilevered  uniform  flat  plate  to  determine  the 
validity  of  the  algorithm,  Henderson  (Ref  3) .  Favorable  results  were 
realized,  as  the  first  eight  (8)  modes  for  eleven  (11)  mass  loading 
configurations  were  predicted  to  within  a  small  margin  of  error  for 


virtually  every  case.  All  but  nine  cases  yielded  errors  of  less  than  10%, 
and  those  other  nine  were  all  less  than  20%.  Whaley  (Ref  14). 

Purpose 

Even  though  this  technique  has  proven  itself  acceptable  in  the  cases 
of  beams  and  flat  plates,  further  validation  is  required  in  the  specific 
area  of  built-up  aircraft  structures.  Since  this  methodology  was 
originally  developed  with  these  types  of  structures  in  mind,  it  is 
necessary  to  determine  whether  this  theory  holds  valid  for  a  structure 
of  much  greater  complexity  than  those  previously  tested,  one  typical  of 
what  can  be  found  on  today's  aircraft.  For  this  purpose,  a  curved,  rib- 
stiffened  panel  from  a  cargo  aircraft  was  selected  for  experimental 
validation.  The  structure  has  two  longitudinal  stiffeners,  and  five 
lateral  stiffeners.  A  complete  modal  survey  was  conducted  to  determine 
the  first  several  natural  frequencies  and  their  respective  mode  shapes. 

The  selected  number  of  modes  to  be  included  in  the  analysis  will  be 
discussed  in  Chapter  II  along  with  the  reasons  for  chosing  it.  The 
modal  analysis  will  be  done  for  the  panel  in  its  unloaded  configuration, 
and  nine  separate  mass  loaded  configurations. 

Following  completion  of  the  above  experimental  testing,  the  unloaded 
modal  data  was  utilized  in  the  algorithm  to  predict  the  natural  frequencies 
for  each  mass  loading  configuration.  The  experimentally  determined  values 
will  then  be  compared  to  the  predicted  ones  to  find  the  percent  error, 
and  ultimately  the  validity  of  the  algorithm.  After  it  was  determined 
that  these  results  were  favorable,  the  system  forced  response  was 
evaluated  for  each  mass  loaded  configuration.  These  will  then  be 
compared  to  the  experimentally  determined  values,  to  ultimately 
discover  the  accuracy  of  this  prediction  method  when  it  is  applied 
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to  complex  structures. 

Previous  Mass  Loading  Techniques 

There  presently  exist  a  number  of  techniques  for  modeling  various  types 
of  structures  with  an  added  mass  load.  Several  of  these  are  discussed 
briefly  by  Whaley  (Ref  14) .  Those  most  closely  related  to  the  technique 
used  in  this  study,  are  that  of  Weissenburger  (Ref  12) ,  and  Pomazal  and 
Snyder  (Ref  6).  Weissenburger  expresses  the  natural  frequencies  of  the 
loaded  structure  in  terms  of  those  of  the  unloaded  structure  which  results 
in  the  need  for  solving  a  transendental  equation  (in  matrix  form)  .  Pomazal 
and  Snyder  basically  employ  Weissenburger 's  technique,  but  utilize  a 
Newton-Raphson  matrix  interaction  technique  to  solve  for  the  modified 
eigenvalues.  The  technique  utilized  in  this  study  differs  in  that  the 
eigenvalue  problem  is  posed  in  terms  of  measured,  unmodified  natural 
frequencies  and  mode  shapes.  The  complexity  of  the  problem  is  directly 
dependent  on  the  desired  number  of  natural  frequencies  and/or  mode  shapes. 
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The  derivation  of  the  algorithm  to  be  used  in  this  study  has  been 


reported  by  Whaley  (Ref  14)  and  is  discussed  in  detail  below.  By 
expressing  the  response  of  a  continuous  structure  in  terms  of  its 
normal  modes,  it  is  possible  to  transform  the  governing  partial 
differential  equation  for  the  structure  into  an  infinite  set  of 
uncoupled  ordinary  different  equations  in  terms  of  the  generalized 
coordinates.  Lagrange's  equation  is: 

d  3T  iV  Q . 

dt  3q^  (1) 


for  a  flat  plate: 


(3) 


where: 

W(x,y,t)  =  I  $1  (x,y)  q±(t)  (4) 

Def ine  6P  =  Virtual  work  of  applied  inertial  loads. 


6P  -  F(x0,y0.t)5W(xo,yo,t)  +  M(x0»y0»t)fiw(x0»y0»t) 


Therefore  the  generalized  force  may  be  expressed  as: 
6P 


«i-  -Mo  j^j(xo»yo>  3E 


<*2qj 


-  Mo  Rx2  M1  <xo»y0)  jii  ^(xo’yo>  ^ 


34-4 


“  “o  %2  if1  ^O’yo5  jii  Sy  -  (xO»yQ>  y1-  (5) 


Combining  Eqs  (1)  through  (5)  yields  an  infinite  set  of  equations, 
the  ith  one  of  which  will  be: 


Mi  h  +  Ki  qi  *  Qi 


(6) 


where  the  following  definitions  hold: 


ba 


Generalized  Mass  =  M  =  ff  ph  <J> . 2  (x,y)  dxdy 

1  oo  1 


(7a) 


Generalized  Stiffness  -  K.  =  //  (ttt - t.  ((^..v*)  +  (*■■»>--) 

-  i  oo  12(l-vz)  3x^  3yz 


32d>4  32*i 


+  2v  a^1  a^T  +  dxdy  {7b) 


Now,  we  can  avoid  the  necessity  of  solving  transcendental  equations 
or  using  an  extremely  complex  finite  element  modelling  technique,  if  we 
can  truncate  our  infinite  set  of  Eqs  (6)  to  a  sufficiently  small  finite 
number.  One  can  accomplish  this  by  assuming  that  a  reliable  result  can 
be  achieved  by  considering  only  the  first  "n"  modes.  Then  Eqs  (6)  can 
be  written  in  matrix  form  as  follows: 

[A]  (q1)  +  [K]  {q}  -  (o)  (8) 
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where  (K)  Is  a  diagonal  matrix  made  up  of  the  generalized  stiffness 
elements  determined  from  Eq  (7b) ,  and  [a]  is  made  up  of  the  following 
elements: 

Asm  -  m  +  Mo  $2m(xo,yo)  +  Mo  r2  (||“  (*o»yo))2  +  Mo  Ry^^^.yo))2 

(8a) 

Amn  -  MoV(xo»yo)*n  <-*o,yo)  +  Mq  r5  (^-n  (xo*y°))  (8b) 

+  41“  (xo’yo>  ^  (x°*yo))  ■  ¥  n 

If  one  now  assumes  simple  harmonic  motion  (i.e.,qBqe  ),  and 
premultiplies  Eq  (8)  by  [k]  \  a  standard  eigenvalue  problem  of 
dimensions  n  by  n  results: 

{4[l]-  [K]'1  [A]}  (i>  =  (0)  (9) 

or 

The  frequency  to  can  be  found  by  setting  the  determinant  of  the  above 
matrix  to  zero. 

Thus,  one  is  faced  with  a  problem  which  requires  a  great  deal 
less  computational  effort  than  if  it  were  necessary  to  derive  a  detailed 
model  for  the  structure.  For  the  complex  structures  that  are  typically 
found  on  aircraft,  the  various  stiffeners,  spars,  ribs,  and  curved 
panels  make  a  continuous  model  out  of  the  question.  Therefore,  the 
most  common  approach  to  this  type  of  structure  has  been  to  formulate 
a  finite  element  model.  The  finite  element  analysis  is,  however,  much 
more  complex,  requires  a  great  deal  more  modelling  and  programming 
effort,  and  utilizes  more  expensive  computer  time  than  solving  Eq  (9) 
for  the  loaded  frequencies. 

In  summary,  to  utilize  this  algorithm,  one  requires  only  knowledge 
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of  the  unloaded  natural  frequencies  and  mode  shapes  (descrete  representa¬ 
tions),  and  the  magnitude,  location  and  radii  of  gyration  of  the  added 


Ill  Experimental  Procedure 

Experimental  testing  was  conducted  to  determine  the  natural 
frequencies  and  discrete  mode  shapes  of  the  unloaded  panel.  These 
data  were  required  by  the  prediction  algorithm  to  approximate  the 
loaded  natural  frequencies.  The  natural  frequencies  of  the  panel 
in  each  of  its  mass  loaded  configurations  also  had  to  be  experimentally 
determined  to  enable  comparison  between  the  actual  and  predicted 
natural  frequencies.  In  this  manner,  the  validity  of  the  algorithm 
for  the  chosen  complex  structure  could  be  determined.  All  testing 
was  conducted  at  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL) . 

All  equipment  that  was  required  is  listed  in  Table  1.  Due  to  problems 
in  acquisition  and  processing  of  data  which  will  be  discussed  below, 
three  separate  test  methods  were  attempted.  The  third  one  of  these 
was  the  only  one  to  yield  acceptable  data.  All  three  of  these  methods 
are  discussed  in  this  chapter  in  some  detail.  This  is  done  with  the 
goal  of  familiarizing  the  reader  with  the  three  methods,  the  aspects 
in  which  each  did  or  did  not  fail,  and  the  reasons  for  their  respective 
failures  or  alternatively  their  successes.  Armed  with  this  knowledge, 
persons  conducting  future  testing  in  this  area  will  be  able  to  avoid 
problems  similar  to  those  encountered  by  the  author  in  the  course  of 
this  study. 

The  panel  was  suspended  in  a  horizontal  position  from  the  ceiling 
by  bungy  cords  attached  to  each  corner.  This  allowed  freedom  of  motion 
of  the  entire  panel  in  any  direction  (termed  a  FREE  -  FREE  mode) .  A 
grid  of  accelerometer  locations  was  then  laid  out  on  the  panel  surface 
as  shown  in  Fig  1.  Three  different  masses  were  chosen  as  mass  loads. 
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TABLE  I  EXPERIMENTAL  APPARATUS 


APPARATUS 

MANUFACTURER 

MODEL 

COMMENTS 

Accelerometer 

Calibrator 

Bruel  and  Kjaer 

B&K  4291 

1  g  peak 

Fourier  Spectrum 
Analyzer 

Hewlett-Packard 

HP  5451B 

25  lb.  Shaker 

Ling  Electronics 

Random  Noise 
Generator 

Hewlett-Packard 

HP  3722A 

flat  spectrum 

Accelerometers  (6) 

Vibrametrics 

Force  Gauge 

Vibrametrics 

208  A03 

1  mV/lb  force 

Amplifiers  (6) 

Intech 

A-2318 

Variable  gain 

Oscilloscope  (2) 

Hewlett-Packard 

HP  1707B 

Universal  Filter 

General  -  Radio 

Bandpass  50-1000  Hz 

Terminal 

Tektronix 

TEK  4014-1 

Copier 

Xerox 

Versatek 

Hard  copies  from 
terminal 

FM  Recorder 

LAR 

LAR  7400 

Voltmeter 

NLS 

LX-2 

Power  Supply 

AFFDL/FBG 

+  15V  DC 

Strobe 

Strobex 

121  A 

Force  Gauge 

Power  Unit 

Piezotronlcs 

480  A 

Low  Frequency 
Oscillator 

Electrodyne 

N  300 
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The  isometric  views  of  these  three  masses  are  shown  in  Figures  2,  3,  and  4 
Then,  three  arbitrary  mass  loading  locations  were  chosen  as  shown  in  Fig  5 
Utilizing  each  mass  in  all  three  locations  yielded  a  total  of  nine  mass 
loaded  configurations.  The  pertinent  data  for  all  nine  of  these  is  given 
in  Table  II.  In  each  of  the  three  separate  testing  procedures  all  of  the 
above-mentioned  conditions  remained  the  same.  Briefly,  the  three  methods 
employed  were: 

1.  Single  point  excitation  using  a  continuous,  random  noise  input 
to  the  25  lb.  shaker.  All  accelerometer  data  was  recorded  on  FM  tape, 
along  with  the  force  input  (from  a  force  gauge) .  The  recorded  data  was 
then  digitized  and  processed  to  yield  Power  Spectral  Density  and  Transfer 
Function  plots. 

2.  The  force  input  to  the  panel  was  the  same  as  above,  but  all 
data  processing  was  done  by  the  HP5451B  Fourier  Spectrum  Analyzer. 

3.  The  input  excitation  was  15  impulse  inputs  at  random  locations. 
Data  processing  was  done  by  the  HP5451B,  and  all  response  data  was  an 
average  Of  the  responses  to  each  of  the  15  Inputs.  These  three  methods  are 
described  in  detail  below. 

Test  Method  1 

The  panel  was  set  up  as  was  previously  discussed.  A  25  lb.  Ling 
Shaker  was  attached  to  one  of  the  two  longitudinal  stiffeners  which  ran 
the  length  of  the  panel,  almost  directly  below  the  first  column  of 
accelerometers,  and  was  slightly  off  center  from  the  middle  lateral 
stiffener.  The  input  signal  to  the  shaker  was  a  continuous,  broadband, 
random  noise  signal.  A  force  gauge  was  placed  between  the  shaker  and  the 
panel.  Due  primarily  to  a  software  problem  which  rendered  the  Hewlett- 
Packard  Fourier  Spectrum  Analyzer  unavailable  when  testing  was  begun,  an 
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Table  II  HASS  LOADING  CONFIGURATION  DATA 


Config¬ 

uration 

Number 

Mass 

Number 

Mass  (Kg) 

Rx  (m) 

m 

Y0  (m) 

1 

1 

0.110 

0.0129 

0.0129 

0.152 

0.279 

2 

1 

0.110 

0.0129 

0.0129 

0.356 

0.356 

3 

1 

0.110 

0.0129 

0.0129 

0.508 

0.533 

4 

2 

0.189 

0.0142 

0.0142 

0.152 

0.279 

5 

2 

0.189 

0.0142 

0.0142 

0.356 

0.356 

6 

2 

0.189 

0.0142 

0.0142 

0.508 

0.533 

7 

3 

0.052 

0.0123 

0.0188 

0.152 

0.279 

8 

3 

0.052 

0.0123 

0.0188 

0.356 

0.356 

9 

3 

0.052 

0.0188 

0.0123 

0.508 

0.533 
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older,  more  conventional  method  was  utilized.  The  response  data  was  collected 
from  five  (5)  accelerometers  (one  column)  at  a  time  (i.e.,  accelerometers 
were  placed  first  on  points  (1,1),  (1,2),  (1,3),  (1,4),  (1,5),  then  moved 
to  (2,1),  (2,2),  (2,3),  (2,4),  (2,5),  and  so  on).  These  points  were  as 
depicted  in  Fig  1  and  the  output  of  the  force  guage  and  all  five  accelerometers 
were  amplified  and  recorded  on  FM  tape.  See  Fig  6  for  a  complete  wiring 
diagram.  Upon  completion  of  thirty  (30)  seconds  of  recorded  data,  the 
accelerometers  were  moved  to  columns  2,  3,  4,  and  5  in  succession,  repeating 
the  above  procedure  each  time.  This  FM  tape  was  then  digitized  and  processed 
to  yield  Power  Spectral  Density  (PSD)  and  Transfer  Function  Plots.  See 
Appendix  E  for  a  sample  of  the  Transfer  Function  Magnitude  and  Phase  Plots 
that  resulted.  From  these  it  was  possible  to  determine  the  frequency  band 
of  interest,  as  it  was  desirable  to  confine  the  analysis  to  approximately 
the  first  ten  (10)*  modes,  which  covered  a  frequency  band  of  50  -  500  Hz 
(the  lower  limit  avoids  some  very  low  frequency  modes  of  the  suspension 
system  which  were  found  by  doing  a  frequency  sine  sweep) .  In  addition  to 
the  Transfer  Function  Plots,  tabular  listings  were  printed  for  both  the 
magnitude  and  phase  of  the  Transfer  Function  values  for  intervals  of  0.5  Hz. 
From  these,  the  natural  frequencies  and  a  discrete  representation  of  each 
mode  shape  was  determined.  This  procedure  was  carried  out  for  the  panel  in 
an  unloaded  configuration  and  each  of  the  nine  mass  loaded  configurations. 

At  this  stage,  a  serious  problem  arose  in  that  the  phase  information 

proved  to  be  erroneous.  Theoretically,  the  phase  angle  between  the  input 

force  and  an  output  acceleration  should  always  be  +90°.  This,  however, 

was  not  the  case.  In  fact,  at  every  frequency  that  had  consistently  yielded 

*  Ten  was  chosen  because  it  is  small  enough  that  the  eigenvalue 
problem  does  not  become  too  unwieldly,  but  as  outlined  in  Whaley 
(Ref  5)  this  is  a  sufficient  number  of  terms  to  provide  a  good 
approximation  of  forced  response. 
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Fig.  6  Wiring  Diagram  -  Test  Method  1 
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a  local  peak  in  the  Transfer  Function  plots,  and  was  therefore  thought  to  be 
a  natural  frequency,  there  were  many  points  at  which  the  phase  angle  was  as 
much  as  40  away  from  +  90  .  Because  of  this  discrepancy,  coherence  function 
plots  were  examined,  which  are  in  essence  a  measure  of  the  reliability  of  the 
data.  The  coherence  function  is  defined  in  Richardson  (Ref  9)  as  follows: 


o 

is  (to)  |2  * 

2 

1  y*  1 

S  (to)s  (to) 

XX  yyV  7 

(10) 

where:  Input  Power  Spectrum  ■  S^ 

(to)  =  X(to)X*(to) 

(11) 

Output  Power  Spectrum  =  S  (to)  =  Y(to)Y*(to) 

(12) 

Cross  Power  Spectrum  =  S 

r  xy 

(to)  =  Y(to)X*(to) 

(13) 

Representing  the  system  in  block  form: 


Fig.  7  Block  Diagram  of  Test  Station 
X(to)  =  Fourier  Transform  of  input  signal 
H(to)  =  System  Transfer  Function 
Z(<o)  =  Fourier  Transform  of  desired  output  signal 

N(to)  =  Fourier  Transform  of  extraneous  noise  (thermal  noise,  60-cycle,  etc.) 
Y(to)  =  Fourier  Transform  of  actual  measured  output  signal 


*  Bar  denotes  algebraic  average  of  the  number  of  measurements 
taken. 
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(14) 


Y(oj)  “  H(td)  +  N(oj) 

S„(b)  =>  Y(to) X*(w)  -  H(s)X(io)  +  N(u>)  X*(w) 
xy 

S  (to)  =  H(co)  S  (to)  +  N(u>)X*(io)  (15) 

xy  xx 

Syy(to)  =  (H(to)X(to)  +  N(to))  (H(to)X(to)  +  N(to))* 

S  (to)  *  |H(<o)  1 2  S  (to)  +  H*(to)  S  (to)  +  H(to)S  (to)  +  S  (to)  (16) 

yy  xx  nx  xn  nn 

S  (to)  =  (H(to)X(to))  (R(co)X(to))* 
zz 

S  (to)  =  |H((o)  j 2  S  (to)  (17) 

XX  XX 

From  Eqs  (10) ,  (15)  and  (16)  it  can  be  seen  that  when  measured  response 
power  is  equal  to  measured  input  power,  y2  =  1,  and  the  noise  power,  N(o)), 
has  no  effect  on  H(to).  However,  when  measured  response  power  is  greater 
than  measured  input  power  because  of  the  contribution  of  some  extraneous 
noise  sources,  y2<  1,  for  all  those  frequencies  where  the  noise  adds  power 
to  the  response  signal.  Therefore,  the  coherence  function  can  be  used  to 
indicate  the  degree  of  noise  contamination  in  Transfer  Function  measurements. 
Since  the  additive  noise  can  be  modelled  as  a  zero  mean  random  process,  its 
effect  on  the  Transfer  Function  results  can  be  reduced  by  averaging  results 
of  several  responses  together  (i.e.,  ensemble  averaging). 

When  the  coherence  function  plots  for  each  grid  point  were  generated,  it 
was  found  that  for  everyone  of  them  in  all  configurations,  the  coherence  was 
dropping  almost  to  zero  at  several  frequencies  indicating  that  the  measure¬ 
ments  had  been  contaminated  by  a  great  deal  of  extraneous  noise  (See  Fig  8) . 
To  alleviate  this  problem,  a  couple  of  alterations  to  parameters  in  the  data 
processing  method  were  tried:  1)  increasing  the  frequency  resolution  and 
therefore  the  transform  time,  and  2)  increasing  the  number  of  transforms  in 

the  averaging  process.  These  both  improved  the  resulting  data,  but  the 
coherence  function  plots  still  indicated  that  a  great  deal  of  unwanted  noise 
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RIB-STIFFENED  PANEL  TESTS  UNLOADED 
(1.3) 


Fig.  8  Sample  Coherence  Function  Plot 
Showing  How  it  Drops  Off  at  Numerous  Places 


r 


was  present  in  the  system. 

The  author  and  some  engineers  at  the  Air  Force  Flight  Dynamics  Laboratory 
(AFFDL)  formulated  the  following  hypothesis  as  to  the  cause  of  the  problem: 
the  unacceptably  large  amount  of  noise  contamination  was  the  result  of  a 
combination  of  two  factors  -  the  amplitude  of  the  input  excitation  force, 
and  the  extremely  low  structural  damping  factor  of  the  test  specimen.  As  a 
result,  an  input  to  the  panel  at,  say,  T  *  0,  would  not  be  damped  out  for  a 
long  period  of  time.  To  understand  fully  how  this  becomes  a  problem,  one 
must  know  how  data  processing  for  continuous  inputs  is  accomplished.  In  the 
case  of  thirty  (30)  second  runs  of  recorded  data,  and  using  the  maximum 
frequency  resolution  (.4242  Hz)  attainable  on  the  processing  equipment  at 
AFFDL,  means  that  Fourier  Transforms  will  be  taken  over  1/.4242  =  2.357 
seconds  intervals  and  twelve  (12)  transforms  can  be  taken.  These  twelve 
transforms  can  then  be  averaged  to  reduce  the  effect  of  extraneous  noise  on 
the  results.  To  yield  the  desired  output,  artificial  damping  is  introduced 
onto  the  recorded  data  in  the  form  of  a  "Hanning  Window".  (See  Fig  9)  In 
this  manner,  one  can  look  at  the  panel  response  over  a  given  "time  window" 
(2.357  seconds)  in  the  30  second  sample,  then  take  the  Fourier  Transforms  of 
the  input  and  the  output  at  a  given  point,  and  generate  the  Power  Spectral 
Density  or  Transfer  Function  plots  for  that  particular  point.  A  similar 
process  must  be  carried  out  for  all  twenty-five  (25)  points.  What  actually 
happens  in  the  case  of  such  a  lightly  damped  structure  is  that  any  transient 
response  which  is  due  to  inputs  prior  to  the  beginning  of  any  given  "time 
window",  enter  that  window  as  noise.  Note  that  this  effect  would  be  absent 
if  one  could  compute  an  "exact"  Fourier  Transform  (over  an  infinite  time 
interval).  The  reliability  of  all  the  data  is  extremely  dependent  upon  the 
Signal- to-Noise  Ratio  (SNR)  being  sufficiently  high  to  yield  good  data. 
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Fig.  9  Hanning  Window  artificially 
damping  response  to  zero 

Therefore,  as  more  and  more  prior  response  enters  each  successive  window  as 
noise,  the  SNR  becomes  worse  and  worse,  and  so  does  the  data.  The  problem 
could  not  be  resolved  using  the  best  possible  values  of  total  data  time, 
frequency  resolution,  and  number  of  transforms  in  the  averaging  process  that 
AFFDL's  data  processing  unit  had  to  offer. 

In  summary,  this  first  test  method  failed  because  the  extremely  low 
structural  damping  factor  of  the  panel  resulted  in  its  responding  to  a  given 
input  for  a  long  period  of  time.  This  combined  with  the  inability  of  the 
data  processing  method  to  allow  sufficient  time  for  the  panel  responses  to 
damp  out  caused  the  SNR  to  grow  worse  with  each  successive  Fourier  Transform 
taken.  This  in  turn  caused  the  data  taken  to  become  worse  and  worse  with  the 
result  that  the  gathered  data  could  not  be  considered  reliable.  Therefore, 
a  new  test  method  was  sought. 

Test  Method  2 

By  this  point  in  time  the  HP5451B,  Fourier  Analyzer,  Hewlett-Packard  (Ref  4), 
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was  ready  for  use.  This  device  utilizes  Fast  Fourier  Transform  and  Ensemble 
Averaging  Techniques  to  calculate  Power  Spectral  Densities  and  Transfer 
Functions  very  rapidly.  The  test  set  up  was  still  exactly  the  same  as  in 
Fig  6,  with  the  exception  that  because  the  Fourier  Analyzer  could  perform 
all  data  processing,  the  FM  Recorder  was  no  longer  required.  The  two  inputs 
to  the  analyzer  were  the  output  of  the  force  gauge  (always) ,  and  the  output 
of  each  one  of  the  twenty-five  accelerometers  in  succession.  The  panel  was 
still  driven  by  the  twenty-five  lb  shaker  excited  by  broadband  random  noise. 

The  HP5451B  was  capable  of  achieving  better  limits  of  the  performance  para¬ 
meters:  longer  transform  time  and  therefore  better  frequency  resolution,  and 
greater  numbers  of  transforms  to  be  used  in  the  averaging  process.  It  was 
hoped  that  this  improved  performance  would  solve  the  problem,  and  yield  some 
reliable  data.  It  was  also  thought  that  although  it  would  not  improve  the 
coherence  at  any  given  point,  the  removal  of  any  bay  modes  from  the  analysis 
would  simplify  the  task  of  forming  the  discrete  mode  shapes.  Bay  modes  are 
modes  which  occur  in  only  one  or  more,  but  not  all,  of  the  bays  created  by 
the  isolating  effect  of  the  stiffeners.  Since  the  author  was  interested  only 
in  global  modes,  or  modes  of  the  entire  panel,  eliminating  bay  modes  from  the 
analysis  altogether  at  the  beginning  of  testing,  would  save  a  great  deal  of 
time  during  the  data  reduction  process.  To  this  end,  a  high  frequency  strobe 
was  obtained.  The  intention  was  to  do  a  sine  sweep  throughout  the  frequency  . 
band  of  interest,  tune  to  each  harmonic,  and  view  the  panel  with  the  strobe 
to  determine  whether  the  mode  was  a  global  one,  or  just  a  mode  of  one  or  more 
bays.  However,  the  amplitude  of  panel  motion  was  very  small  and  therefore  it 
was  extremely  difficult  to  view  the  panel  modes.  The  strobe  was,  as  a  result, 
ineffective.  The  Fourier  Analyzer  was  also  incapable  of  resolving  the  problems 
of  unreliable  data.  This  was  once  again  determined  by  examining  the  coherence 
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function  plots.  The  Fourier  Analyzer  was  capable  of  computing  these  using 
Eq  (10) ,  and  then  plotting  them.  These  revealed  the  same  problem  as  was 
experienced  earlier.  The  noise  power  at  several  frequencies  was  so  high 
that  the  coherence  dropped  off  very  nearly  to  zero  at  some,  and  to  as  low 
as  0.5  at  many  others.  Therefore  it  was  concluded  that  the  data  was  not 
reliable  and  could  not  be  used.  As  a  result,  a  third  experimental  procedure 
was  required. 

Test  Method  3 

This  third  method,  in  which  the  method  of  data  analysis  was  developed 
by  Brown  (Ref  7)  proved  to  be  effective,  and  was  the  one  finally  used. 

First,  a  location  for  a  "reference  accelerometer"  was  chosen.  This  location 
remained  the  same  throughout  all  testing.  It  was  discovered  that  an 
accelerometer  location  near  the  center  of  the  panel  yielded  better  results 
than  a  location  close  to  an  edge.  The  amplified  output  from  this  reference 
accelerometer  was  one  of  the  inputs  to  the  Fourier  Analyzer.  The  other 
input  was  the  amplified  output  of  an  accelerometer  at  one  of  the  twenty-five 
data  points  on  the  grid  (Fig  1) .  A  complete  wiring  diagram  of  this  experi¬ 
mental  set  up  is  shown  in  Fig  10. 

The  panel  was  then  excited  by  tapping  (impulse  input)  the  surface  once  at 
each  of  fifteen  (15)  arbitrary  locations.  After  each  individual  impulse, 
sufficient  time  was  allowed  for  the  panel  response  to  dampen  out  (an 
indicator  light  on  the  Spectrum  Analyzer  notified  the  tester  when  the 
response  had  decreased  to  a  negligible  level) .  For  each  excitation  the 
Fourier  Analyzer  was  used  to  compute  the  Auto  Power  Spectrum  of  the  response 
at  the  reference  point,  the  Auto  Power  Spectrum  of  the  response  at  a  given 
point,  and  the  Cross  Power  Spectrum  between  the  two.  A  stable  average  of 
each  was  then  computed  over  the  fifteen  inputs.  The  computed  parameters  are 
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defined  as  follows: 


Output  Auto  Power  Spectrum  of  Reference  =  Srr(u)  =  R(w)R*(u>) 
Output  Auto  Power  Spectrum  of  Data  Point  ■  Syy(u>)  -  Y(to)Y*(o)) 

Cross  Power  Spectrum  *  Syr(o>)  a  R(ai)Y*(u)) 

_  , ,  .  .  .  *n  ^n-1 

Stable  Average  ■  A  ■  A  ,  +  - 

n  n-1  n 

Then,  at  the  end  of  the  fifteen  excitations  a  Modal  Assurance 
function  (Brown) ,  was  computed  which  is  defined  as  follows: 


(18) 

(19) 

(20) 

(21) 

Criteria  (MAC) 


MAC 


is  (t*0 
I  yr 


2* 


Srr(“>Syy(w)  (22) 

The  difference  between  this  function,  and  the  coherence  function  defined  in 
Eq  (10)  is  that  both  auto  spectrums  are  of  responses  at  two  points  (there  is 
no  power  spectrum  of  the  force  input) ,  and  the  cross  power  spectrum  is  that 
between  two  output  responses  instead  of  between  a  force  input  and  an  ouput 
excitation  as  in  Eq  (10).  At  all  natural  frequencies,  the  magnitude  and 
phase  relationships  between  the  responses  at  any  two  points  on  the  panel 
surface  will  be  the  same,  as  for  a  given  natural  frequency  there  is  a 
corresponding  mode  shape.  Therefore,  at  these  modal  frequencies,  the 
magnitude  of  the  Cross  Power  Spectrum  squared  will  be  equal  to  the  product 
of  the  individual  Auto  Spectra.  It  follows  that  at  these  frequencies  the 
value  of  the  MAC  function  will  be  1.0.  Now,  at  all  non-modal  frequencies 
the  relationship  between  the  responses  at  any  two  points  on  the  surface  of 
the  panel,  is  a  zero  mean  random  process  (Brown).  Therefore,  for  each  of 
the  fifteen  excitations,  the  relationship  between  the  Reference  and  a  given 
point  will  be  different,  and  the  average  of  the  Cross  Power  Spectrum  over 
those  fifteen  will  be  small  compared  to  1.  Thus,  the  value  of  the  MAC 
function  at  all  frequencies  that  are  not  natural  frequencies  will  be  small. 

*  Bar  (-)  over  values  in  Eq  (22)  denotes  stable  average  of  those  parameters. 


(See  Fig  11  for  a  sample  MAC  function) .  Also  computed  at  the  completion  of 
the  fifteen  impulse  inputs,  was  the  Transfer  Function.  In  the  case  of  this 
type  of  testing,  this  Transfer  Function  was  not  defined  in  the  sense  of  a 
classical  system  Transfer  Function  (which  is  output  response  divided  by 
input  response),  but  was  defined  as  follows: 

H(<u)  - 

Srr  (a))  (23) 

This  is  a  measure  of  the  relative  response  at  a  given  point  as  compared 
with  the  response  at  the  Reference  point. 

In  accordance  with  the  programming  for  this  overall  process,  once  all 
computations  had  been  completed,  the  Fourier  Analyzer  displayed  the  MAC 
Function  plot.  By  viewing  this,  the  user  could  determine  whether  or  not 
the  data  taken  was  acceptable.  If  the  MAC  function  was  a  "clean"  one  as 
shown  in  Figure  11,  i.e.  MAC  =  1.0  at  the  natural  frequencies  and  zero 
elsewhere,  the  data  was  good.  Whereas,  if  there  were  many  values  between 
zero  and  1.0,  and  the  value  remained  at  0.0  very  seldom,  the  data  was 
considered  erroneous.  When  it  was  decided  that  the  data  was  good,  the 
user  gave  the  "CONTINUE"  command,  and  the  Power  Spectral  Density  (PSD) 
plot  (a  measure  of  the  relative  power  in  each  frequency  throughout  the 
band  of  interest) ,  and  the  tabulated  values  of  the  MAC  and  Transfer  Functions 
were  outputted  to  the  TEK  4014-1  scope  in  succession,  and  sent  to  the  Versatek 
hard  copier.  See  Figures  12,  13,  and  14  for  samples  of  each.  This  process 
had  to  be  repeated  for  all  twenty-five  accelerometer  locations,  for  each  of 
the  ten  separate  configurations.  The  PSD  plots  for  each  data  point  of  the 
unloaded  configuration  are  included  in  Appendix  A.  Using  this  procedure 
the  author  was  able  to  generate  the  required  two  hundred  fifty  sets  of 
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Fig.  14  Transfer  Function  (Real  and  Imaginary) 
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data,  all  were  determined  to  be  reliable  on  the  basis  of  the  MAC  function. 

It  should  be  noted  here  that  acquiring  reliable  data  required  re-testing  at 
some  points  where  the  MAC  function  indicated  that  the  data  was  not  good. 

This  bad  data  could  have  been  caused  by  several  factors;  allowing  insufficient 
time  for  panel  response  to  damp  out  between  impulse  excitations,  causing  a 
system  overload  by  striking  the  panel  too  hard  on  one  or  more  of  the  15 
excitations  during  one  test  run,  or  double-tapping  -  mistakenly  striking 
the  panel  twice  on  a  single  impulse  input.  From  this,  it  is  necessary  to 
realize  that  extreme  care  must  be  exercised  when  taking  data,  if  one  expects 
even  this  method  to  yield  reliable  results. 

In  summary,  after  determining  that  data  processing  limitations  precluded 
the  use  of  a  broad-band  random  force  input,  a  method  using  15  discrete  impulse 
inputs  was  used.  An  HP5451B  Fourier  Analyzer,  Hewlett-Packard  (Ref  4)  was 
used  to  perform  all  data  processing.  Using  a  reference  accelerometer  which 
remained  in  a  constant  location  and  an  accelerometer  at  one  of  the  25  data 
points  as  the  inputs  to  the  Fourier  Analyzer,  the  individual  Power  Spectra 
and  the  Transfer  Function  were  computed.  From  these,  the  natural  frequencies 
and  a  discrete  representation  of  the  mode  shapes  were  determined  utilizing 
the  procedure  outlined  in  the  section  that  follows. 
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IV  Computational  Procedure 


ta  Reduction 

This  section  contains  a  description  of  the  procedure  that  was  employed 
extract  the  required  natural  frequencies  and  mode  shapes  from  all  of  the 
a  that  was  collected.  To  yield  this  needed  information  one  must  pick  out 
;  location  of  modes  in  the  frequency  domain,  and  then  find  the  relative 
:nitudes  and  phase  of  the  displacement  of  each  point,  at  each  mode.  This 
t  be  done  for  all  ten  configurations.  By  examining  the  tabulated  values 
ihe  MAC  function  for  a  particular  point  and  for  a  given  configuration,  it 
possible  to  locate  the  natural  frequencies  with  relative  ease,  as  they 
indicated  by  the  occurrence  of  several  consecutive  values  of  the  MAC 
ng  equal  to  or  very  close  to  one  (see  Fig  13,  as  an  example  of  how  the  MAC 
etion  indicates  a  mode  check  frequencies  close  to  168Hz  and  It  can  be 
n  that  there  are  6  consecutive  values  all  approximately  *  1) .  Then  by 
.king  the  corresponding  frequencies  at  the  tabulated  Transfer  Function 
ies,  it  was  possible  to  pinpoint  the  natural  frequency.  This  was  done 
.amining  the  listed  values  in  the  vicinity  of  a  chosen  frequency,  which 
given  in  terms  of  real  and  imaginary  components.  By  realizing  that  the 
e  angle  between  the  two  acceleration  responses  (reference  and  test 
it)  at  a  mode  must  be  either  0°  or  180°,  it  was  necessary  only  to  find 
frequency  where  the  imaginary  component  was  zero  or  closest  to  zero, 
e  this  was  done,  the  relative  position  of  that  point,  in  the  positive  or 
tive  direction  was  given  by  the  real  component  of  the  transfer  function 
e  at  the  determined  frequency.  This  procedure  was  carried  out  for  all 
lty-five  points  for  each  configuration.  This  yielded  the  natural 
uencies  and  a  discrete  representation  of  the  corresponding  mode  shapes 
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for  the  first  nine  (9)  modes  (all  the  modes  that  existed  between  0  and  500  Hz) 
for  all  10  configurations. 

An  example  of  the  discrete  mode  shape  plots  is  shown  in  Fig  15,  and  the 
plots  for  the  first  nine  modes  of  the  unloaded  configurations  are  included  in 
Appendix  8. 

Description  of  Computer  Programs 
General 

Since  both  Eq  (9)  and  the  computation  of  the  generalized  elements  required 
derivatives  of  the  mode  shapes,  continuous  representations  of  the  discrete 
mode  shapes  previously  formed,  must  somehow  be  generated.  This  was  accomplished 
using  a  bi-cubic  spline  curve  fitting  approach.  The  bi-cubic  spline  is  a 
mathematical  algorithm  which  fits  a  cubic  polynomial  between  sets  of  discrete 
data  points,  and  yields  continuous  first  and  second-order  derivatives  at  these 
data  points.  Also  required  for  the  computation  of  the  Generalized  Masses  were 
the  thickness  and  mass  density  of  the  panel.  Since  one  of  the  objectives  of 
this  prediction  technique  was  to  enable  analysis  of  a  complex  structure  in  a 
relatively  simple  manner,  without  having  to  utilize  a  complex  model  for  the 
structure  (i.e.,  finite  element  approach),  the  specimen  was  modelled  as  a 
uniform  panel  for  computational  purposes.  A  geometric  averaging  approach 
was  used  to  take  the  stiffeners  into  account.  The  surface  area  of  the 
stiffeners  was  averaged  over  the  surface  area  of  the  panel  to  yield  a  con¬ 
stant  thickness  which  was  somewhat  greater  than  that  of  the  panel  surface 
alone.  The  mass  of  the  stiffeners  was  also  taken  into  account  by  then 
calculating  an  effective  panel  mass  density  using  the  known  overall  panel 
mass,  surface  area,  and  the  averaged  thickness. 

Once  this  was  done,  the  generalized  mass  and  stiffness  were  computed 
for  each  mode  of  each  configuration,  using  Program  1  described  below.  These 
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Fig.  15  Discrete  Mode  Shape  Plot 


These  values  were  then  used  In  Eq  (9) ,  to  solve  for  the  predicted  natural 
frequencies  (of  each  mass  loaded  configuration),  using  Program  2  described 
below. 

Program  1 

This  program  had  been  largely  developed  previously,  but  was  adapted  and 
improved  for  the  purpose  of  this  study.  A  complete  listing  of  this  computer 
program  may  be  found  in  Appendix  D.  Utilizing  only  unloaded  natural  frequency 
and  discrete  mode  shape  data,  this  program  made  use  of  existing  International 
Mathematical  and  Statistical  Libraries  (IMSL)  subroutines  to  compute  a  bicubic 
spline  approximation  to  curve-fit  a  continuous  mode  shape  to  the  discrete  data. 
The  user  controls  the  accuracy  of  this  approximation  by  chosing  the  number  of 
points  in  both  the  X  and  Y  direction  at  which  he  wishes  the  values  to  be 
computed.  In  the  case  of  this  study,  the  author  chose  to  divide  the  total 
length  in  the  X  and  Y  directions  into  forty  (40)  parts.  This  means  that  the 
value  of  each  mode  shape  was  computed  at  a  total  of  1600  points.  Obviously, 
the  finer  one  chooses  the  mesh  size  for  these  computations,  the  more  accurate 
will  be  the  approximation  of  the  continuous  mode  shape,  but  the  computer  time 
and  therefore  expense  also  increases  drastically.  Therefore,  a  certain  amount 
of  engineering  judgement  was  required  in  the  determination  of  the  desired 
mesh  size.  Once  the  continuous  mode  shape  approximations  had  been  computed, 
they  were  used  in  Eq  (7a) ,  and  two  separate  integral  equations  were  solved, 
one  for  the  surface  of  the  specimen  itself,  and  one  for  the  discrete  mass 
load.  The  two  results  were  then  summed  to  yield  the  Generalized  Mass  values 
for  each  of  the  first  nine  modes,  for  each  one  of  the  nine  mass  loading 
configurations.  The  Generalized  Stiffness  values  were  then  computed  using 
the  following  formula: 
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(24) 


Vi 
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Where  the  to^'s  are  the  unloaded  natural  frequencies.  The  "j"  subscript 
represents  the  configuration  number.  Note  that  Eq  (7b)  was  not  used 
to  compute  the  Generalized  Stiffness  values  due  to  the  large  errors  that 
result  from  squaring  the  second  derivative  of  the  approximated  mode  shapes. 

A  complete  listing  of  the  Generalized  Mass  and  Stiffness  values  computed 
by  this  program  are  included  in  Appendix  D. 

Program  2 

Again,  much  of  this  program  had  been  previously  written  by  Whaley  (Ref  3) 
for  the  test  conducted  on  the  flat  plate.  However,  it  was  necessary  to 
adapt  and  improve  it  due  to  the  greater  complexity  and  scope  of  this  study. 

This  program  solves  Eq  (9)  for  the  predicted  mass-loaded  natural  frequencies. 

As  can  be  seen  from  Eqs  (8a)  and  (8b) ,  aside  from  the  nine  (9)  Generalized 
Mass  and  Stiffness  Values  for  each  configuration  (from  Program  1)  other  data 
values  were  required  as  inputs  to  this  program.  These  were,  the  value  of 
the  added  mass  (Kg),  the  location  of  the  added  mass,  (Xo,Yo),  as  shown  in 
Fig  5,  and  the  radii  of  gyration  of  the  added  mass,  (Rx,Ry),  See  Table  II 
for  all  this  data.  A  sample  calculation  of  the  radii  of  gyration  for  one  of 
the  masses  loads  used  is  shown  in  Appendix  E.  A  complete  listing  of  this 
program  is  included  in  Appendix  D  and  the  outputs  showing  the  predicted 
loaded  natural  frequencies  (as  percentages  of  the  unloaded  ones)  are  included 
at  Appendix  D. 

Program  3 

This  program  is  really  Program  1  with  a  couple  of  deletions,  and  a  large 
addition  written  by  the  author  which  allows  the  user  to  output  three-dimensional 
plots  of  as  many  of  the  continuous  mode  shapes  as  is  desired.  This  was  used 
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largely  to  clear  up  any  ambiguities  that  arose  when  deciding  which  of  the 
experimentally  determined  mode  shapes  for  all  of  the  mass  loading  configurations 
matched  up  with  those  from  the  unloaded  case.  Due  to  some  modes  being  very 
close  to  one  another  and  some  being  shifted  upwards  in  frequency,  it  was 
necessary  to  plot  the  mode  shapes  and  compare  them  to  the  unloaded  ones. 

In  this  manner,  it  was  possible  to  determine  the  correspondence  between  the 
loaded  modes  shapes  and  the  unloaded  ones,  where  any  doubt  had  previously 
existed.  As  one  example  of  this,  a  mode  existed  at  238.22  Hz  in  configuration  6. 
This  is  at  a  slightly  higher  frequency  than  the  nearest  mode  in  the  unloaded  case 
(237. 74  Hz).  The  mode  shapes  are  almost  exactly  the  same  leaving  no  question 
that  they  were  the  same  mode.  A  second  example,  is  that  a  mode  was  found  at 
257.5  Hz  in  configuration  7,  and  here  again,  matching  up  the  mode  shapes 
confirmed  this  mode  was  indeed  the  same  as  the  one  at  252.10  Hz  in  the 
unloaded  case.  In  each  case,  two  views  of  the  same  mode  shape  were  generated. 
Examples  of  these  plots  are  shown  in  Figures  16  and  17.  A  complete  listing 
of  the  program  used  to  make  these  plots  is  included  in  Appendix  D.  The  two 
views  plotted  of  the  first  nine  mode  shapes  for  the  unloaded  configuration 
plus  those  from  the  two  examples  mentioned  above,  are  included  in  Appendix  C. 
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Fig.  16  View  1  -  4th  Mode  -  Unloaded  Configuration 


Fig.  17  View  2  -  4th  Mode  -  Unloaded  Configuration 


A 


V  Results 


A  summary  of  the  test  results  and  significant  findings  is  presented  in 
this  section.  Table  III  gives  all  the  experimentally  determined  natural 
frequencies  for  the  unloaded  and  loaded  configurations.  Some  very  interesting 
results  were  immediately  noticeable.  The  first  mode  at  167.93  Hz  in  the 
unloaded  case  was  suppressed  (did  not  show  up  at  all)  in  the  case  of 
configurations  1,  2,  3,  and  4  which,  because  it  is  the  first  mode,  should 
have  a  significant  impact  on  the  forced  response,  and  on  the  error  incurred 
in  the  forced  response  predictions.  This  was  confirmed  by  plotting  the  first 
two  or  three  mode  shapes  in  each  case,  and  matching  them  up  with  those  for 
the  unloaded  case.  Also,  since  there  was  a  chance,  albeit  remote,  that  the 
referenced  accelerometer  had  been  placed  on  a  node  of  the  first  mode  for  each 
configuration,  which  would  cause  the  mode  to  seem  to  disappear,  the  panel  was 
re-tested  in  these  four  configurations  with  a  new  reference  location.  However, 
the  results  were  unchanged.  Also  interesting  was  the  fact  that  for  some 
configurations,  some  modes  were  actually  shifted  upward  somewhat  infrequency. 
Modes  5  and  6  of  configuration  7  are  examples  of  this.  Again,  these  results 
were  confirmed  by  plotting  these  particular  mode  shapes  and  comparing  them  to 
the  corresponding  unloaded  shapes.  In  each  case,  there  was  no  question  that 
the  modes  matched  up  as  they  are  shown  in  Table  III.  This  second  phenomenon 
seems  counter  intuitive  but  a  possible  explanation  for  it  is  offered  here. 

It  has  been  found,  Talmadge  (Ref  8),  that  in  the  case  of  structures  such  as 
the  panel  used  in  this  study  that  for  a  given  global  mode,  the  individual 
bays  (formed  by  the  stiffeners)  often  do  not  all  vibrate  at  exactly  the  same 
frequency  or  amplitude.  One  of  the  bay  frequencies  could  be  predominant  in 
the  unloaded  case,  but  the  location  of  the  mass  load  makes  one  of  the  other 
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Table  III 


NATURAL  FREQUENCIES  -  EXPERIMENTAL 


167.93 

180.50 

191.22 

207.31 

237.74 

252.12 

173.83 

184.23 

201.05 

237.52 

253.09 

174.08 

187.62 

204.00 

225.31 

243.32 

171.40 

181.42 

210.80 

240.70 

251.08 

173.90 

183.10 

211.50 

237.11 

252.53 

167.70 

179.53 

187.14 

202.58 

223.92 

246.32 

169.71 

180.72 

199.30 

208.27 

238.22 

252.20 

143.10 

175.09 

184.83 

204.30 

238.22 

257.50 

153.70 

177.02 

189.63 

205.21 

230.94 

245.82 

162.02 

177.43 

182.71 

205.64 

237.22 

251.43 

r 


bays  the  dominant  one.  Since  it  is  not  unusual  for  the  bay  frequencies  to 
be  4  or  5  Hz  apart  it  is  very  possible  to  find  a  loaded  mode  that  is  a  few 
Hz  higher  in  frequency  than  the  corresponding  unloaded  one.  The  phenomenon 
of  modal  suppression  however,  could  not  be  theoretically  explained  by 
the  author  or  any  of  the  test  engineers  at  AFFDL,  although  it  had  been 
encountered  in  previous  testing  there.  In  spite  of  this,  due  to  two 
factors,  the  data  was  considered  to  be  correct.  (1)  Since  there  had  been 
so  little  previous  testing  in  this  area  with  a  structure  as  complex,  and 
therefore  so  little  knowledge,  exactly  what  to  expect  is  an  unknown 
parameter.  (2)  Since  such  care  had  been  taken  in  gathering  the  data, 
and  the  indications  from  the  MAC  functions  were  that  the  data  was  good, 
the  author  was  confident  that  the  data  was  reliable. 

Tables  IV  through  XII  show  the  comparison  between  the  experimentally 
determined  and  the  predicted  natural  frequencies,  along  with  the  resulting 
per  cent  error  of  each  predicted  value.  This  is  shown  separately  for 
each  configuration.  Excluding  the  four  occurrences  of  complete  modal 
suppression,  in  all  but  one  case  the  per  cent  error  was  well  below  10% 
and  the  one  case  where  this  is  not  so,  the  error  is  well  below  20%.  Note 
that  in  Table  VI  (configuration  3)  and  in  Table  IX  (configuration  6)  the 
predicted  natural  frequencies  for  modes  7  and  8  are  reversed  in  order. 

That  is,  the  predictions  reflect  that  the  8th  mode  is  actually  shifted 
below  the  7th  mode  in  frequency.  This  is  not  what  happens  in  reality  and 
therefore,  the  prediction  of  the  eighth  natural  frequency  in  both  of  these 
cases  is  considered  a  major  failure  on  the  part  of  the  algorithm,  even 
though  the  per  cent  errors  in  these  cases  are  still  only  5.50%  and  7.21% 
respectively.  Another  interesting  phenomenon  occurred  that  should  be 
noted  here:  that  was  the  occurrence  of  modes  in  some  mass  loaded 
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Table 

IV  DATA  RESULTS  - 

CONFIGURATION  1 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

(  ) 

167.67 

(  ) 

2 

173.83 

179.16 

3.07 

3 

184.23 

190.94 

3.64 

4 

201.05 

206.88 

2.90 

5 

237.52 

237.09 

-  0.18 

6 

253.09 

246.89 

-  2.45 

7 

279.79 

276.23 

-  1.27 

8 

292.77 

292.14 

-  0.22 

9 

364.89 

364.18 

-  0.19 

Table  V 

DATA  RESULTS  - 

CONFIGURATION  2 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

(  ) 

167.79 

(  ) 

2 

174.08 

180.15 

3.49 

3 

187.62 

191.13 

1.87 

4 

204.00 

206.43 

1.19 

5 

225.31 

237.37 

5.35 

6 

243.32 

251.86 

3.51 

7 

263.80 

279.61 

5.99 

8 

282.45 

291.95 

3.36 

9 

363.08 

362.76 

-  0.09 

Table 

VI  DATA  RESULTS  - 

CONFIGURATION  3 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

(  ) 

166.46 

(  ) 

2 

171.40 

176.36 

2.89 

3 

181.42 

191.25 

5.42 

4 

210.80 

206.10 

-  2.23 

5 

240.70 

235.81 

-  0.69 

6 

251.08 

251.51 

0.17 

7 

280.00 

278.95 

-  0.38 

8 

292.22 

276.16 

-  5.50 

9 

340.95 

344.80 

1.13 
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Table  VII  DATA  RESULTS  -  CONFIGURATION  4 


MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

(  ) 

167.69 

(  ) 

2 

173.90 

179.03 

2.95 

3 

183.10 

190.83 

4.22 

4 

211.50 

206.69 

-  2.27 

5 

237.11 

236.83 

-  0.11 

6 

252.53 

245.04 

-  2.97 

7 

279.63 

274.65 

-  1.78 

8 

291.88 

292.03 

0.05 

9 

341.60 

363.76 

6.49 

Table 

VIII  DATA  RESULTS  - 

CONFIGURATION  5 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

167.70 

167.79 

0.054 

2 

179.53 

180.10 

0.56 

3 

187.14 

190.01 

1.53 

4 

202.58 

207.20 

2.31 

5 

223.92 

237.32 

5.98 

6 

246.32 

251.79 

2.22 

7 

273.65 

279.44 

2.12 

8 

291.19 

291.77 

0.20 

9 

362.81 

361.34 

-  0.41 

Table 

IX  DATA  RESULTS  - 

CONFIGURATION  6 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

169.71 

166.22 

-  2.06 

2 

180.72 

175.34 

-  2.98 

3 

199.30 

191.23 

-  4.05 

4 

208.27 

205.63 

-  1.27 

5 

238.22 

235.20 

-  1.27 

6 

252.20 

251.18 

-  0.40 

7 

279.95 

278.44 

-  0.54 

8 

292.20 

271.14 

-  7.21 

9 

360.20 

338.70 

-  6.11 
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Table  X  DATA  RESULTS  -  CONFIGURATION  7 


MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

143.10 

167.77 

17.24 

2 

175.09 

180.20 

2.91 

3 

184.83 

191.16 

3.44 

4 

204.30 

206.98 

1.31 

5 

238.22 

237.33 

-  0.37 

6 

257.50 

250.45 

-  2.73 

7 

279.88 

278.70 

-  0.42 

8 

292.76 

292.48 

-  0.096 

9 

361.71 

364.39 

0.74 

Table  XI 

DATA  RESULTS  - 

CONFIGURATION  8 

PERCENT  ERROR 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

1 

153.70 

167.80 

9.17 

2 

177.02 

180.41 

1.93 

3 

189.63 

191.23 

0.86 

4 

205.21 

206.86 

0.81 

5 

230.94 

237.50 

2.85 

6 

245.82 

251.99 

2.52 

7 

273.73 

279.88 

2.26 

8 

281.02 

292.36 

4.04 

9 

362.84 

363.69 

0.25 

Table 

XII  DATA  RESULTS  - 

CONFIGURATION  9 

MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

162.02 

166.73 

2.92 

2 

177.43 

177.77 

0.21 

3 

182.71 

191.27 

4.69 

4 

205.64 

206.53 

0.45 

5 

237.22 

236.48 

0.30 

6 

251.43 

251.51 

0.044 

7 

279.51 

279.43 

-  0.025 

8 

292.29 

281.42 

-  3.72 

9 

361.48 

351.59 

-  2.74 
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configurations,  that  did  not  show  up  at  all  in  the  unloaded  case.  For 
example,  there  was  a  mode  at  202.7  Hz  in  configuration  3  and  one  at  201.0  Hz 
in  configuration  4.  However,  the  author  offers  the  following  possible 
explanation  of  this  phenomenon:  the  "extra"  mode  would  have  been  there 
in  the  unloaded  case,  but  was  so  close  (in  frequency)  to  another  mode  that 
they  appeared  as  one.  But  these  two  particular  mass  loading  configurations 
simply  caused  them  to  be  shifted  apart  sufficiently  to  make  one  distinguishable 
from  the  other.  This  phenomenon  is  discussed  in  Richardson  (Ref  9)  and  it 
was  determined  that  this  is  a  frequent  occurrence  in  testing  of  this  kind, 
in  a  private  conversation  with  the  chief  test  engineer  at  the  Air  Force 
Flight  Dynamics  Laboratory  Talmadge  (Ref  8) .  Therefore,  the  fact  that 
there  was  no  predicted  counterpart  for  these  "extra"  modes  was  not 
considered  an  error  on  the  part  of  the  prediction  algorithm. 

In  summary,  the  algorithm  predicted  the  mass-loaded  natural  frequencies 
to  within  a  very  small  amount  of  error,  excluding  those  unexplainable 
occurrences  of  modal  suppression.  So,  with  the  exception  of  those  four 
cases,  the  results  were  quite  favorable. 
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VI  Discussion 


General 

In  74  of  81  cases  the  natural  frequencies  of  the  panel  were  predicted 
with  less  than  10%  error,  and  one  other  was  less  than  20%  error.  Of  the 
remaining  six  cases,  four  were  the  occurrences  of  complete  modal  suppression 
which  the  algorithm  has  no  way  of  predicting.  This  would  become  an  obvious 
shortcoming  in  the  event  that  this  phenomenon  occurred  on  a  large  number 
of  occasions  for  one  particular  structure.  But  based  on  the  fact  that  in 
the  case  of  this  particular  structure  it  only  happened  a  total  of  four 
times  and  never  occurred  more  than  once  for  a  given  mass  loaded  config¬ 
uration,  the  event  of  numerous  such  occurrences  is  unlikely.  The  last 
two  cases  were  those  of  the  modal  "flip-flopping"  in  configurations  3  and 
6,  discussed  in  the  previous  chapter.  These  errors  were  more  than  likely 
due  to  three  factors: 

1.  For  the  purposes  of  the  prediction  program,  the  panel  was  modelled 
as  a  flat  plate  of  uniform  thickness  and  density,  which  was  the  result  of 

•;  the  geometric  smearing  of  the  stiffeners  discussed  earlier. 

2.  The  bi-cubic  approximation  of  the  continuous  mode  shapes  were 
computed  only  over  the  surface  area  covered  by  the  grid,  and  not  the 
entire  surface  area  of  the  panel.  This  was  due  to  the  fact  that  only 
some  and  not  all  of  the  subroutines  employed  in  this  process  were  capable 
of  approximating  values  beyond  the  boundaries  of  the  gathered  data. 

3.  Addition  of  the  point  mass  loads  caused  warping  of  the  mode 
shapes.  The  loaded  mode  shapes  are  not  exactly  the  same  as  the  unloaded 
ones,  and  since  the  unloaded  mode  shapes  are  used  in  the  prediction  cal¬ 
culations,  some  error  naturally  arises. 

i 
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The  first  nine  natural  frequencies  were  predicted  under  nine-  different 


mass  loading  conditions.  Therefore,  a  broad  array  of  loading  conditions 
was  examined.  There  does,  however,  remain  some  doubt  in  the  author's  mind 
that  this  spectrum  of  mass  loading  conditions  was  broad  enough  to  conclude 
with  any  certainty  that  this  method  is  effective  for  all  arbitrary  mass 
loading  configurations.  A  careful  examination  of  Table  III  reveals  that 
the  largest  change  in  any  of  the  loaded  natural  frequencies  was  only  18.7%. 
(Coincidentally,  this  was  the  case  with  the  largest  error  between  measured 
and  predicted  which  was  equal  to  17.24%.)  Most  of  the  loaded  natural 
frequencies  were  altered  by  only  5%  or  less  from  their  unloaded  counter¬ 
parts.  In  light  of  this,  the  error  margins  shown  in  Tables  IV  through 
XII  may  not  be  quite  as  impressive  as  they  look  at  first  glance.  It  is 
suspected  that  if  the  unloaded  frequencies  were  changed  by  a  significantly 
greater  amount,  that  these  error  margins  may  worsen  in  direct  proportion. 
The  fact  that  the  loaded  frequencies  were  changed  by  the  small  amount 
mentioned  is  directly  attributable  to  the  size  of  the  mass  loads  used. 

In  defense  of  this,  fairly  small  mass  loads  were  chosen  for  two  reasons: 

1.  The  panel  was  suspended  from  the  ceiling  by  bungy  cords  in  a 
configuration  somewhere  between  FREE  -  FREE  and  SIMPLY  SUPPORTED.  Due 
to  this,  it  was  thought  that  it  would  not  require  large  mass  loads  to 
alter  the  vibration  characteristics  by  an  easily  measurable  amount. 

2.  On  the  recommendation  of  the  chief  test  engineer  at  AFFDL,  who 
thought  that  even  the  largest  mass  used  in  this  study  might  be  too  large, 
the  heavier  masses  that  the  author  originally  planned  to  use  were  removed 
from  consideration  during  this  study.  In  an  effort  to  prove  or  disprove 
the  hypothesis  that  the  error  was  directly  related  to  the  amount  the 
natural  frequencies  were  modified,  the  author  (at  the  time  of  this  writing) 
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returned  to  the  lab  and  reset-up  all  the  test  equipment  as  shown  in  Fig  10. 
Two  new  masses  of  significantly  larger  size  than  the  others  were  used  in 
one  test  run  each.  The  first,  a  .259  Kg  rectangular  hexahedron,  was  placed 
in  location  1,  denoted  configuration  10,  as  shown  in  Fig  5.  The  second, 
a  .592  Kg  solid  cylinder,  was  placed  in  location  3,  denoted  configuration  11. 
Each  was  tested  in  the  same  manner  as  that  outlined  in  Chapter  III  under 
"Test  Method  3".  Then,  all  the  necessary  computations  were  performed  on 
the  computer  to  yield  the  predicted  natural  frequencies  and  the  gathered 
data  was  reduced  to  yield  the  actual  values.  The  results  of  these,  along 
with  the  per  cent  error  incurred  are  included  in  Tables  XIII  and  XIV. 

In  order  to  determine  whether  or  not  this  additional  testing  accomplished 
what  it  was  hoped,  these  new  results  were  compared  to  those  found  using 
the  smaller  masses  in  the  same  locations.  The  results  from  configuration  10 
were  compared  to  those  from  configurations  1,  4,  and  7  (all  masses  in 
location  1)  and  the  results  from  configuration  11  were  compared  to  those 
from  configurations  3,  6,  and  9  (all  masses  in  location  3).  Unfortunately 
even  the  heavier  masses  used  did  not  modify  the  natural  frequencies  as  much 
as  it  was  hoped.  As  a  result,  this  new  information  was  not  sufficient  to 
enable  any  concrete,  indisputable  conclusions  to  be  drawn.  However,  some 
valuable  information  was  gleened  from  it.  Comparing  the  configurations  in 
the  manner  stated  above,  it  was  seen  that  for  both  of  these  two  new  cases, 
the  error  margins  are  generally  higher  than  any  of  the  three  other  respective 
cases  for  each.  Accordingly,  most  of  the  natural  frequencies  were  altered 
somewhat  more  than  before.  This  indicated  that  the  accuracy  of  the  algorithm 
diminishes  hand  in  hand  with  increasing  magnitude  of  the  mass  load  with 
respect  to  panel  mass.  This  performance  degradation  is  not  a  desirable 

property,  as  It  means  that  the  algorithm  could  only  be  used  in  cases  where 
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Table  XIII  DATA  RESULTS  -  CONFIGURATION  10 


MODE 

ACTUAL  FREQUENCY 

PREDICTED  FREQUENCY 

PERCENT  ERROR 

1 

(  ) 

167.59 

(  ) 

2 

173.13 

177.72 

2.65 

3 

183.17 

189.83 

3.64 

4 

201.09 

205.62 

2.25 

5 

239.25 

235.45 

-  1.59 

6 

256.05 

240.50 

-  6.07 

7 

280.42 

270.96 

-  3.37 

8 

293.71 

291.91 

-  0.61 

9 

332.18 

362.07 

9.00 

Table  XIV  DATA  RESULTS  -  CONFIGURATION  11 


MODE  ACTUAL  FREQUENCY  PREDICTED  FREQUENCY  PERCENT  ERROR 


1 

2 

3 

4 

5 

6 

7 

8 
9 


(  ) 
170.25 
180.18 

208.42 

239.42 
252.90 
280.23 
293.11 
311.41 


162.03 

167.20 

190.86 

200.81 

229.67 

241.59 

273.83 

236.81 

295.48 


(  ) 

-  1.79 
5.92 

-  3.65 

-  4.07 

-  4.47 

-  2.28 
-19.21 
-  5.12 
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the  mass  loads  are  relatively  small,  such  that  they  cause  little  change  in 
the  natural  frequencies.  This  area  needs  to  be  examined  in  relation  to 
the  forced  response,  to  discover  how  much  error  will  result  in  calculating 
the  predicted  forced  response  for  a  given  amount  of  error  in  the  prediction 
of  the  natural  frequencies.  Such  an  examination  would  yield  a  yardstick  for 
determining  how  much  error  in  the  natural  frequencies  is  acceptable.  This 
is  discussed  briefly  in  the  section  that  follows. 

Forced  Response 

This  section  develops  a  method  of  calculating  the  Power  Spectrum  of  the 
forced  response  for  a  hypothetical  case  where  the  forcing  function  is  a 
white  noise  process.  Then,  since  comparing  two  power  spectra  to  determine 
any  error  incurred  would  require  plotting  in  4  -  dimensions,  the  RMS  (Root 
Mean  Square)  value  of  the  forced  response  is  calculated,  giving  a  good 
indicator  of  the  overall  difference  between  two  forced  responses.  A  "best 
case"  and  "worst  case"  error  model  is  then  examined  for  a  constant  error 
in  all  predicted  natural  frequencies.  It  should  be  noted  here  that  no 
comparison  of  the  actual  versus  predicted  forced  response  could  be  done 
during  the  course  of  this  study  due  to  the  manner  in  which  the  experimental 
testing  was  conducted.  The  force  input  excitation  finally  used  after 
encountering  the  problems  outlined  in  Chapter  III  of  this  study  was  fifteen 
impulse  excitations  of  random  strength  and  location.  Since  these  inputs 
were  not  one  of  the  parameters  measured  (response  at  two  points,  reference 
and  one  data  point,  were  measured)  there  was  no  known  way  to  calculate  the 
actual  system  transfer  function  or  to  compute  statistical  information  (such 
as  the  Powpr  Spectrum)  of  the  input.  Without  this  information,  the 
experimental  forced  response  could  not  be  computed  and  therefore  could 
not  be  compared  to  any  predicted  forced  response. 
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d.  Integrate  over  all  frequencies  to  yield  the  RMS  value  of  the 
responses  (also  called  average  power) : 

E{x(t)Y(t) }  -  R  (o) 

xy 
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Error  Model 

a.  Worst  case.  The  worst  error  was  in  configuration  11,  and  it  was 
-19.212.  So,  if  it  is  assumed  that  all  predicted  frequencies  are  20%  low: 
let  n=*9  as  in  this  study 

9  Mx*y) 

R  (o)-l.SvHN  ,Z.— J — 2—2 - J! - 

xy  oo  j=l  (>8  uj2 


9  <Mxn»y  >Mx>y> 

1.5ir  M  N  ,Z.  — J — 2_2 — L 


0  0  64u  2 

J 
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Development .  From  Bogdanoff  (Ref  1)  we  have: 

S  (x.y.w)  -  I  {[M>(x,y)^i(x,y)dxdy  S.(a))ij)  (x,y)]|H(a»)  |2>  (25) 

j“l  oo  J  1  J 

where:  P(x,y)  °  panel  loading  as  a  function  of  x  and  y 
<f^(x,y)  *  unloaded  mode  shapes 

S^(m)  =  Power  Spectrum  the  Forced  Input 

H(m)  =»  System  Transfer  Function 
n  »  number  of  modes  being  considered 

S  (x,y,w)  =  Power  Spectrum  of  the  response  as  a  function  of  x,  y 
w  frequency 

a.  Substituting  in  expression  for  the  System  Transfer  Function: 

1 


H(u>) 


“j  2-u>2+2  j  Cmcoj 


where  =  panel  natural  frequencies 

5  ■  damping  factor  chosen  =*  .05 

ba 

n  ^^x»yH.(x.y)dxdyS  (u>)«t>  (x,y) 

S  (x,y,w)  -  - L— - i - i- - }  (26a) 

^  1  “  “2)2  +  (2(*05)“wj)2 

b.  Substitute  in  expression  for  Panel  Loading  and  Integrate: 

for  point  mass  loading: 

P(x»y)  =  M S(x  ,y  ) 
o  o  o 

where  M  ■  magnitude  of  mass  load 
o 

6  =*  dirac  delta  function 

«v(*,y.U)  -  j  {Mb) 

(o)2-U2)2  +  (Orl  OJOJj)2 

c.  Modelling  random  noise  input  as  a  white  noise  process: 

S  («#)  -  N 

1  o 


,  and 
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1  9 

{1.5M  N  l 
.64  °  o  ]»1 


<j>.(x  ,y  H  (x,y) 

_i _ o  Is  2i _ 


.  64 


R  (o) 
xy  correct 


=  1.563  R  (o) 

xy  correct 

Therefore,  if  the  natural  frequencies  were  all  predicted  20%  below  their 
true  values,  the  predicted  RMS  value  of  the  response  is  56.3%  higher  than 
the  true  value  (note  that  this  discounts  any  error  that  is  a  result  of 
mode  shape  warping  due  to  the  mass  load  which  will  result  because  unloaded 
mode  shapes  are  used  in  this  prediction) .  This  amount  of  error  is 
obviously  unacceptable. 

b.  Best  case.  Again,  looking  at  configuration  11,  the  smallest  error 
incurred  was  -1.79%.  Therefore,  assuming  that  the  predictions  are  all 
only  2%  low: 


R  (o) 
xy 


9 

1.5*  MN  E, 
o  o  j=l 


*.(x0,y0H.(x,y) 

(.98^)2 


9 

-  1.5n  M  N  ,1. 

o  o  j=l 


<Mx  ,y  H.(x,y) 

-10  0-1 

.96u)j2 


1 

.96 


9 
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o  o  j=] 


♦  .<*  »y  ><t>4(x,y) 
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This  time  only  a  4.2%  error  results  which  for  most  applications  would  be 
an  acceptable  value. 

c.  General  Result 

E  =  {[ - - -  -  1]  x  100}  %  (28) 

(1+e)  2 

E  =  overall  error  in  the  prediction  of  RMS  value  of  response 
e  =  %  error  of  natural  frequency  predictions  expressed  as  a  decimal 
Note  that  this  hold  only  if  all  predictions  are  in  error  by  the  same 
amount  and  all  are  either  high  ojr  low.  If  any  of  these  conditions  do 
not  hold,  then  Eq  (24)  must  be  solved  implicitly  to  determine  the  exact 
error  (requires  a  computer  program,  as  the  calculation  involves  operations 
on  3-dimensional  mode  shapes) .  It  is  also  important  to  notice  that  in 
the  most  of  the  11  configurations  finally  tested,  some  of  the  errors  are 
positive,  and  some  are  negative.  In  cases  like  that,  some  large  errors 
may  be  acceptable  as  they  will  have  a  cancelling  effect  if  the  signs  are 
opposite  from  one  another. 

Using  this  "yardstick"  to  examine  the  error  in  the  output  forced 
response  of  the  nine  original  configurations,  the  two  worst  cases  would 
be  6  (all  errors  negative  -  I.e.,  all  predictions  low)  and  8  (all  errors 
positive) . 

i.  Configuration  6.  The  highest  errors  were  -7.21,  and  -6.11,  but 
these  occurred  on  the  two  highest  order  modes  where  the  magnitude  of  the 
natural  frequencies  lessens  their  importance  in  the  summation  (24) . 
Therefore,'  consider  the  average  error  =  -2.88%  (simply  an  algebraic 
average  of  all  the  errors  in  configuration  6). 


E  =,  [ - 1 - -  _  i]  x  (100) 

(l-.0288)^ 

=*  6.011  % 

which  is  an  acceptable  error. 

ii.  Configuration  8.  Here  the  largest  error  =  9.17 %  and  is  on  the 
first  or  most  important  term  so  to  yield  an  indication  of  what  to  expect 
a  weighted  average  is  computed  by  doubling  this  worst  error.  Now  average 
error  =  +4.17% 

E  =  [ - i - ,  -  1]  x  (100) 

(1  +  .0417) 

-  -7.846% 

which  is  also  acceptable. 

One  must  remember  that  these  are  only  good  indications  of  what  the  actual 
error  would  be.  Since  all  the  errors  are  differenct,  Eq  (27)  must  be 
solved  implicitly  to  find  exact  error.  However,  these  figures  provide 
a  reliable  figure  upon  which  some  conclusions  can  be  made.  These  are 
discussed  in  the  section  that  follows.  It  is  also  important  not  to  forget 
that  these  calculations  determine  only  the  error  resulting  from  the  inexact 
predictions  of  the  natural  frequencies,  and  do  not  consider  those  due  to 
mode  shape  warping.  This  error  can  only  be  found  by  comparing  the  pre¬ 
dicted  forced  response,  and  the  experimentally  determined  one.  As  was 
discussed  earlier,  this  was  not  possible  during  the  course  of  this  study, 
and  must  be  left  to  future  research.  Also,  from  an  examination  of  Eq  (27) 
it  must  be  noted  that  since  the  first  mode  is  divided  by  the  lowest  number 
(lowest  frequency)  its  suppression  is  likely  to  have  significant  effect 


on  the  forced  response  which  will  cause  a  large  error  in  the  prediction 
as  the  algorithm  does  not  predict  its  disappearance.  Exactly  how  large 
this  error  is,  must  also  be  left  to  future  studies. 
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VII  Conclusions  &  Recommendations 


Conclusions 

The  following  conclusions  can  be  stated  regarding  the  use  of  the 
algorithm  as  a  means  of  predicting  mass  loaded  natural  frequencies  of 
complex  structures: 

1)  The  algorithm  is  both  effective  and  accurate  for  predicting  the 
mass  loaded  natural  frequencies  when  the  lumped  masses  are  of  Hie  same 
relative  magnitudes  as  those  used  during  the  course  of  this  study  (up 
to  20%  of  total  panel  mass) . 

2)  On  the  basis  of  the  forced  response  error  model  calculations  the 
algorithm  also  yields  acceptably  small  errors  in  estimating  the  forced 
response. 

3)  Some  mass  loading  configurations  cause  the  suppression  of  the 
first  mode  which  will  have  a  significant  effect  both  on  the  forced  response 
of  the  structure  and  on  the  error  of  the  prediction. 

4)  Results  indicate  that  as  the  relative  magnitudes  of  the  mass 
loads  increase  so  do  the  errors  in  predicting  the  natural  frequencies 
and  therefore  the  forced  response.  So,  the  applications  of  the  algorithm 
are  limited  to  comparatively  small  loads. 

5)  The  rate  at  which  these  errors  increase  seems  to  depend  not  only 
upon  the  magnitude  of  the  added  mass,  but  also  on  its  location  on  the  panel 
surface  and  radii  of  gyration.  This  is  borne  out  by  the  fact  that  in  the 
case  of  the  rectangular  259  Kg  mass  errors  increased  in  almost  exactly 
direct  proportion  to  the  increase  in  mass,  but  for  the  cylindrical 

.592  Kg  mass,  the  average  error  increased  at  only  slightly  greater  than 
half  the  rate  of  increasing  mass. 
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Recommendations 


It  is  recommended  that  further  studies  be  conducted  to  determine  the 
following: 

1)  An  exact  measure  of  the  accuracy  of  this  method  as  a  means  of 
predicting  the  forced  response  of  complex  structures,  and  the  effect  of 
suppressed  modes.  To  accomplish  this,  a  different  test  method  than  the 
one  finally  employed  here  should  be  used.  There  are  two  recommended  test 
procedures  outlined  in  Appendix  F.  Both  are  capable  of  supplying  all 
information  necessary  to  compute  the  experimental  force  response. 

Briefly,  these  methods  are: 

a.  Use  exactly  the  same  test  procedure  as  was  outlined  earlier 
under  Test  Method  1,  but  with  a  new  data  processing  method  (called  "Overlap 
Processing")  that  was  not  available  at  the  time  of  this  study. 

b.  A  classic  approach  using  an  impulse  hammer  with  a  force  gauge 
to  measure  the  input  force,  and  an  accelerometer  to  measure  output  response. 
The  Fourier  analyzer  does  all  data  processing. 

Once  the  experimental  response  has  been  computed,  the  predicted  force 
response  can  be  calculated  using  the  development  given  in  the  previous 
Chapter,  and  then  the  two  can  be  compared.  In  these  future  tests,  it  is 
also  suggested  the  grid  layout  of  accelerometer  locations  be  extended  to 
the  edges  of  the  test  specimen.  This  will  eliminate  the  error  incurred  by 
excluding  the  edges  of  the  panel  from  the  bi-cubic  spline  continuous  mode 
shape  calculations,  as  discussed  in  the  previous  chapter. 

2)  The  limitations,  in  terms  of  the  relative  magnitudes  of  the  mass 
loads  on  the  capability  of  the  algorithm  to  predict  the  natural  frequencies 
to  within  acceptable  errors.  In  other  words,  since  indications  are  that 


the  accuracy  of  the  algorithm  decreases  as  the  size  of  the  mass  loads  is 
increased,  there  is  a  limit  to  the  algorithm's  effectiveness.  Beyond  this 
limit,  a  different  method  will  have  to  be  used. 


Bibliography 


1.  Bogdanoff,  J.L.  and  Goldberg,  J.E.,  "On  the  Bernoulli  -  Euler  Beam 
Theory  with  Random  Excitation",  Journal  of  the  Aero/Space  Sciencies, 

May  1969,  pp.  371-376.  '  " 

2.  Eringen,  A.C.,  "Response  of  Beams  and  Plates  to  Random  Loads",  Journal 
of  Applied  Mechanics,  March  1957,  pp.  46-51. 

3.  Henderson,  D.A. ,  Whaley,  P.W. ,  and  Brown,  D.L. ,  "Test  Report  on  Modal 
Survey  of  a  Massloaded,  Cantilevered  Flat  Plate",  AFFDL/FBG/78-3, 

April  1978. 

4.  Hewlett  Packard,  "HP5451  B  Fourier  Analyzer",  Operating  and  Service 
Manuals,  Vol.  1-8,  Hewlett  Packard  Co.,  January  1974. 

5.  Lee,  J.  and  Whaley,  P.W.,  "Prediction  of  the  Angular  Vibration  of 
Aircraft  Structures",  Journal  of  Sound  and  Vibration,  Vol.  49(4), 

1976,  pp.  541-549. 

6.  Pomazal,  R.J.,  and  Snyder,  V.W.,  "Local  Modifications  of  Damped  Linear 
Systems",  AIAA  Journal,  Vol.  9,  No.  11,  November  1971,  pp.  2216-2221. 

7.  Private  conversation.  Prof.  D.  Brown,  Head  Mechanical  Eng.  Dept., 
University  of  Cincinnati,  14  July  1979. 

8.  Private  conversation,  Mr.  R.  Talmadge,  Test  Engineer,  AFFDL,  16  August 
1979. 

9.  Richardson,  M. ,  "Modal  Analysis  Using  Digital  Test  Systems",  Seminar 

on  Understanding  Digital  Control  and  Analysis  in  Vibration  Test  Systems, 
Published  by  the  Shock  &  Vibration  Information  Center,  1976. 

10.  Thomson,  W.T.,  Vibration  Theory  and  Applications,  Englewood  Cliffs, 

New  Jersey,  1965. 

11.  Timoshenko,  S.,  Yound,  D.H.,  and  Weaver,  Jr.,  W. ,  Vibration  Problems 
in  Engineering,  4th  Edition,  New  York,  1974. 

12.  Weissenburger,  J.T. ,  "Effect  of  Local  Modifications  on  the  Vibration 
Characteristics  of  Linear  Systems",  Journal  of  Applied  Mechanics, 

Vol;  35,  June  1968,  pp.  327-332. 

13.  Whaley,  P.W.,  "Calculation  of  Natural  Frequencies  and  Mode  Shapes  of 
Mass  Loaded  Aircraft  Structures",  Shock  and  Vibration  Bulletin,  No.  48, 
Pt,  3,  1978,  pp.  13-20. 


62 


14.  Whaley,  P.W. ,  "Prediction  of  the  Change  in  Natural  Frequency  of  a 
Cantilevered  Flat  Plate  with  added  Lumped  Mass",  Not  Yet  Published. 

15.  ,  "Overlap  Processing",  Spectral  Dynamics  Corporation, 
DSP-022,  June  1978. 


63 


Fig.  A-l  PSD  Plot  -  Point  (1,1) 
Fig.  A-2  PSD  Plot  -  Point  (1,2) 
Fig.  A-3  PSD  Plot  -  Point  (1,3) 
Fig.  A-4  PSD  Plot  -  Point  (1,4) 
Fig.  A- 5  PSD  Plot  -  Point  (1,5) 
Fig.  A-6  PSD  Plot  -  Point  (2,1) 
Fig.  A- 7  PSD  Plot  -  Point  (2,2) 
Fig.  A-8  PSD  Plot  -  Point  (2,3) 
Fig.  A-9  PSD  Plot  -  Point  (2,4) 
Fig.  A-10  PSD  Plot  -  Point  (2,5) 
Fig.  A-ll  PSD  Plot  -  Point  (3,1) 
Fig.  A-12  PSD  Plot  -  Point  (3,2) 
Fig.  A-13  PSD  Plot  -  Point  (3,3) 
Fig.  A-14  PSD  Plot  -  Point  (3,4) 
Fig.  A-15  PSD  Plot  -  Point  (3,5) 
Fig.  A-16  PSD  Plot  -  Point  (4,1) 
Fig.  A-l 7  PSD  Plot  -  Point  (4,2) 
Fig.  A- 18  PSD  Plot  -  Point  (4,3) 
Fig.  A- 19  PSD  Plot  -  Point  (4,4) 
Fig.  A-20  PSD  Plot  -  Point  (4,5) 
Fig.  A-21  PSD  Plot  -  Point  (5,1) 
Fig.  A-22  PSD  Plot  -  Point  (5,2) 
Fig.  A-2 3  PSD  Plot  -  Point  (5,3) 
Fig.  A-24  PSD  Plot  -  Point  (5,4) 
Fig.  A-2 5  PSD  Plot  -  Point  (5,5) 


1000  8000  3000  4000  :  5000 


* 

r> 

i 


i 


9  <9 

VT»  VO 

I  I 


I 


PSD  Plot  -  Point  (1,4) 


68 


e8- 


l®  " 


1000  5090  3900  4000  :  5900 


*A0-A079  856  AIR  FORCE  INST  OF  TECH  *RI 5HT-PATTERS0N  AFB  OH  SCHOO— ETC  F/S  20/4 

THE  PREDICTION  OF  HASS  LOADED  NATURAL  FREQUENCIES  ANO  FORCED  RE--ETC(U) 
OEC  79  L  8  SLENESA 

UNCLASSIFIED  AFIT/SAE/AA/79D-5  NL 


Appendix  B 

Discrete  Mode  Shape  Plots  -  Unloaded  Case 
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Fig.  B-8  Discrete  Mode  Shape  -  Mode  8 


Fig.  C-l  (a)  View  1  -  Mode  1  -  Unloaded 

(b)  View  2  -  Mode  1  -  Unloaded 

Fig,  C-2  (a)  View  1  -  Mode  2  -  Unloaded 

(b)  View  2  -  Mode  2  -  Unloaded 

Fig.  C-3  (a)  View  1  -  Mode  3  -  Unloaded 

(b)  View  2  -  Mode  3  -  Unloaded 

Fig.  C-4  (a)  View  1  -  Mode  4  -  Unloaded 

(b)  View  2  -  Mode  4  -  Unloaded 

Fig.  C-5  (a)  View  1  -  Mode  5  -  Unloaded 

(b)  View  2  -  Mode  5  -  Unloaded 

Fig.  C-6  (a)  View  1  -  Mode  6  -  Unloaded 

(b)  View  2  -  Mode  6  -  Unloaded 

Fig.  C-7  (a)  View  1  -  Mode  7  -  Unloaded 

(b)  View  2  -  Mode  7  -  Unloaded 

Fig.  C-8  (a)  View  1  -  Mode  8  -  Unloaded 

(b)  View  1  -  Mode  8  -  Unloaded 

Fig.  C-9  (a)  View  1  -  Mode  9  -  Unloaded 

(b)  View  2  -  Mode  9  -  Unloaded 

Fig.  C-10(a)  View  1  -  Mode  at  238.22  Hz  -  Configuration  6 

(b)  View  2  -  Mode  at  238.22  Hz 


Fig.  C-ll(a)  View  1  -  Mode  at  257.2  Hz  -  Configuration  7 
(b)  View  2  -  Mode  at  257.2  Hz 
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Fig.  C-7  (a)  View  1  -  Mode  7;  (b)  View  2  -  Mode  7  (280.02  Hz) 
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Fig.  C-9  (a)  View  1  -  Mode  9;  (b)  View  2  -  Mode  9  (365.00  Hz) 
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Fig.  C-10  (a)  View  1  -  Mode  at  238.22  Hr  -  Configuration  6 
(b)  View  2  -  Mode  at  238.22  Hr  -  Configuration  6 
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Computer  Program  Listings 
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»?r»f.RAH  si  Mb  ?*/?*  ••r»l 
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c 
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fc 


THiS  PI..1GRAH  CAl**Hir*  1MI  CiMr»ILt4ttf  MASSIf-*.)  AN  0  lift  114*1  |  ?|0 

ff  irt  nc  ?vi6S»  or  r»CM  Mont  f  i,e  a;  »a%mt  cchfigu.  ai  ion*  *s  the  user 

1MIIAS  OMA. 

AiH  W  AM  310*01  !'lr*lT  *MU*  PUT*!  APE  *»  t'HUt*  1  APff  *0'ITM»Y  ,Plf IlE>0> 

PIi-.n  sjuii  r,Mli»  »r.'n»  ,r  spi  #1  ■ 

NUM1£<  nr  MAT*'hAI  rk!OU*Nr.lfS  4riM6  ttAMWH), 

UiHU‘5  I'M  CUSM  T»2»l  t  VIST  ,TtM  »SI4*  »l  »CStl 1M<S,W,41 

roLPtno.  niM>TMr  'wnri  tjr  r*io  p?ihts  ai  vhich  oma  mas  ta<em. 

JtltiVI  *•  r.»ro  “01UIS  Ml  t  AHn  T  111*  CltOHS  RESPECTIVELY* 

Dlhll.SlO'4  l)«PHlsoi-.ifAt»|MU'CI?2bM 

GX(lt)  H<,Tf«»or  Points  t«  *  4NJ  Y  OtRtCUOMS  AT  MMiCM  VALUES 

or  Tl’ l  liciwin  STLIMf  ARE  OLSi^CP. 

IfOV.UOX  VfCTJR  Jr  Lf  MrtMIlMT  ♦4t*MXAMf  »|  *VA*t5  M*-4,VNY-f>  ,rOR 

iu‘:m>uu  ie  oajTtu. 
oihtsion  wnt?rfl»|ptr»'  »*P0st6l 

vcctot  OF  L  ENC  Ml  »•  HA  *  (•’.»  MY  I  *MAt 1 1 Nt»1 1 *3* ‘  J»MTl  #POR 
Summutihe  leoiou. 

fa  lit  riy  NY  HATRI>  O'  fNCTlOM  VA.UCS* 
tlMtl.SlOM  WOK<(  2i  wf> 

WoUfHO**  VECrot  or  trUGlMIMAXt  INT*II*  1*llY»*rOR  SNIROUfltft 
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COI’MtN/PARM/E,f  .  •  • 
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COMMl  N/Mtl0C5/20SB,t  *Y#fSPLIN 
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DATA  ri,II,lY#m,NM.i 
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Ill!  il'fli  »■  ><'.'•5'.,.  211,.  312,. 121. -.2' 3/ 

Sill  HI5!!**1’ ’I’1  t.t.n, M.I/ 
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1111  ISS^IS* 11  ,,rs"  *">  '-•«»».  1’4, -.14*  >.,24, -.145/ 

HI*  *«*«*»'•  "14*0,. ,21, -.4,4/ 

0414  ICOSPIS.il, 1*21, 2M/. 311,. *|?,..,44,,2|”-.l»2/ 
0414  (I0SPI4.II  ,l.|,«ll-.in,.l.,||.1|t|„|!j„lj„ 
0411  (COS  PIG,  1 1  ,1  ,u,  llt/-.|M,,1|I,,,|l|,.,|I||..,M/ 

Sill  KX*!!J,!,,!*U*1'  l/*,tOl2*»C2>.10l*,-.C7t.,-.2»*( 
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0414  lf.«l2PI»,tl,l*l„l3)2..J32„JJ3,.,H,,..|,f,.r„/ 
Sill  IClljP  17,11,1*11,1'  I  /-.24l>,.l2t,-il«,,-.,2  j,  .,.4/ 
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0*11  llOSP14.ll  ,1 -IS,?- I (.123, t»i,..,3, ...  to, 

Dili  IC0jP14.il, 1*21, 2  |2.2t„, l»i, 1*22, (!.;,.  .5J/ 
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Fig.  E-l  Sample  Transfer  Function  Plot  -  Method  1 
Fig.  E-2  Sample  Phase  Plot  -  Methodl 

Fig.  E-3  Sample  Radii  of  Gyration  Calculation 


(3.5) 
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Fig.  E-l  Sample  Transfer  Function  Plot  -  Method  1 
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Appendix  F 

Recommended  Testing  Procedures 

There  are  two  test  methods  that  would  lend  themselves  well  to  a  study 
seeking  to  determine  both  the  accuracy  of  the  algorithm  in  predicting 
loaded  natural  frequencies,  and  the  mass  loaded  forced  response.  The  choice 
of  methods  is  largely  the  user's,  but  will  also  be  dependent  upon  the 
capabilities  of  the  data  processing  unit,  as  the  second  method  covered 
here  requires  knowledge  of  a  very  new  procedure. 

Procedure  1 

This  one  is  a  very  classical  approach,  but  supplies  the  exact 
information  required,  and  is  simple  to  set  up  and  use.  The  test  specimen 
could  be  configured  as  desired  and  as  fine  a  grid  as  is  judged  necessary 
must  be  laid  out  on  the  panel  surface.  Then,  with  accelerometers  in  the 
proper  locations,  the  panel  is  excited  by  striking  it  with  a  hammer  at 
a  single  point.  The  inputs  to  the  HP  5451B  (which  is  capable  of  performing 
all  necessary  data  processing) ,  are  the  amplified  output  of  one  of  the 
accelerometers,  and  the  amplified  output  of  a  force  gauge  which  is  attached 
to  the  hammer.  This  must  be  carried  out  for  each  point  at  which  data  is 
desired.  The  fourier  analyzer  would  have  computed  and  stored  the  input 
force  power  spectrum,  the  output  response  power  spectrum  (from  which  can  be 
found  natural  frequency  and  mode  shape  data) ,  and  the  system  transfer 
function  for  each  point.  This  is  all  the  data  required  for  computation 
of  the  forced  response. 
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Procedure  2 


This  method  is  just  a  repeat  of  Test  Method  1  described  in  Chapter  III. 
The  set  up  can  be  exactly  the  same  and  the  input  force  can  again  be  a 
broadband  random  force  input.  What  allows  one  to  now  use  such  a  method  and 
avoid  the  problems  described  earlier,  is  a  new  data  processing  procedure 
called  Overlap  Processing,  (Ref.  15).  It  removes  the  limitation  on  total 
recorded  data  time,  and  on  the  length  of  the  time  window  over  which  each 
transform  is  taken.  Each  frequency  response  function  is  computed  by 
taking  many  overlapping  transforms  and  averaging  them  all  to  yield  the 
desired  result.  In  such  a  manner  sufficient  time  may  be  allowed  that  the 
response  signal  to  a  given  excitation  can  decay  away  to  zero.  Also, 
sufficient  numbers  of  the  input  and  output  spectrum  can  always  be  computed 
to  allow  adequate  reduction  of  system  statistical  fluctuation,  through 
ensemble  averaging.  It  should  be  noted  here  that  this  technique  is  new, 
that  AFFDL  did  not  learn  of  it  until  approximately  two  weeks  prior  to  this 
writing,  and  it  still  had  not  been  implemented  at  the  time  of  this  writing. 
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generalized  mass  and  generalized  stiffness  for  each  mass  loaded  configuration 
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Block  20  (cont) 

were  computed,  and  the  natural  frequencies  for  each  configuration  were 
computationally  predicted  using  the  aforementioned  algorithm.  The  theoretical 
and  experimental  results  were  then  compared  to  determine  the  amount  of  error 
incurred  in  the  approximation  technique.  The  theory  of  how  to  ultimately 
determine  the  overall  forced  response  of  the  specimen  was  discussed  and  an 
error  model  was  developed  to  enable  an  examination  of  the  reliability  of 
the  algorithm  in  predicting  forced  response.  Recommendations  concerning 
future  test  procedures,  areas  requiring  further  study,  and  the  use  of  the 
algorithm  were  made. 
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